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Abstract: 

ln lH~[U the first author published an algorithm for the conversion of analytic func- 

oc 

tions for which derivative rules are given into their representing power series a,kZ k 

k=0 

at the origin and vice versa, implementations of which exist in Mathematica Jl9[ , 
(s. §), Maple @ (s. g) and Reduce @ (s. @). 

One main part of this procedure is an algorithm to derive a homogeneous linear 
differential equation with polynomial coefficients for the given function. We call this 
type of ordinary differential equations simple. 

Whereas the opposite question to find functions satisfying given differential equations 
is studied in great detail, our question to find differential equations that are satisfied 
by given functions seems to be rarely posed. 

In this paper we consider the family F of functions satisfying a simple differential 
equation generated by the rational, the algebraic, and certain transcendental func- 
tions. It turns out that F forms a linear space of transcendental functions. Further 
F is closed under multiplication and under the composition with rational functions 
and rational powers. These results had been published by Stanley who had proved 
them by theoretical algebraic considerations. 

In contrast our treatment is purely algorithmically oriented. We present algorithms 
that generate simple differential equation for f + g, f -g, f or (r rational), and / ox p ^ q 
(p, q £ JNo), given simple differential equations for /, and <?, and give a priori estimates 
for the order of the resulting differential equations. We show that all order estimates 
are sharp. 

After finishing this article we realized that in independent work Salvy and Zimmer- 
mann published similar algorithms. Our treatment gives a detailed description of 
those algorithms and their validity. 
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1 Simple functions 

Many mathematical functions satisfy a homogeneous linear differential equation with polynomial 
coefficients. We call such an ordinary differential equation simple. Also a function / that satisfies 
a simple differential equation is called simple. The least order of such a simple differential equation 
fulfilled by / is called the order of /. 

Technically there is no difference between a differential equation with polynomial coefficients, 
and one with rational coefficients, as we can multiply such a differential equation by its common 
denominator, so we may consider the case that the coefficients are members of the field K[x] where 
K[x] is one of Q[x], R[x], or C[x]. 

Examples of simple functions are 

1. all rational functions p/q (p, q polynomials) are simple of order 1: pq f + (pq 1 — qp') / = 0, 

2. all algebraic functions / are simple (see e. g. |2[-|3[,[f8|, and [ 10 1 ) of the order of /, 

3. the power function x a satisfies the simple differential equation x f' — a f = 0, 

4. the exponential function e x satisfies the simple differential equation /' — / = 0, 

5. the logarithm function lnx satisfies the simple differential equation x f" + /' = 0, 

6. the sine function sinx and the cosine function cosx satisfy the simple differential equation 

r+/ = o, 

7. the inverse sine function arcsinx and the inverse cosine function arccosx satisfy the simple 
differential equation (x 2 — 1) /" + x f' = 0, 

8. the inverse tangent function arctan x and the inverse cotangent functions arccot x satisfy the 
simple differential equation (1 + x 2 ) f" + 2 x f = 0, 

9. the inverse secant function arcsecx and the inverse cosecant function arccscx satisfy the 
simple differential equation (x 3 — x) f" + (2x 2 — 1) /' = 0, 

10. the error function erf x satisfies the simple differential equation f" + 2 x f = 0. 

Note that for rational functions with rational coefficients, algebraic functions given by a rational 
coefficient equation, the power function for a £ Q, and in all other examples we have K[x] = Q[x]. 



§ 10.4), 
0], 0), 



Also many special functions are simple with K[x] = Q[x]: Airy functions (see e. g. [R 
Bessel functions (see e. g. [Q], Ch. 9-11), all kinds of orthogonal polynomials (see e. g. 
and functions of hypergeometric type (see || ) . 

Not all elementary transcendental functions, however, are simple, the simplest example of which 
probably is the tangent function tanx, compare |l5|, Example 2.5. 



Theorem 1 The functions tanx and secx do not satisfy simple differential equations. 

Proof: It is easily seen that the tangent function f(x) = tanx satisfies the nonlinear differential 
equation 

f = l + f 2 - (1) 
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Differentiation of (jl]) yields after further substitution of (jlj) 

r = (i+/ 2 ) / = 2// / = 2/(i+/ 2 ), 

and inductively we get representations of = Pk{f) (k G N) by polynomial expressions Pk in /. 
Assume now, the tangent function would satisfy the simple differential equation 

n 

P k f [k) = (P* polynomials) , (2) 

fc=0 

then we can substitute by arriving at an algebraic identity for the tangent function. This 
obviously is a contradiction to the fact that the tangent function is transcendental, and therefore 
a differential equation of type (j2|) cannot be satisfied. 

Let us now consider f(x) = secx. In this case we may establish the nonlinear differential 
equation 

f' 2 = f~f 

from which, by differentiation, the second order differential equation 

f" = 2/ 3 - / (3) 
follows. If now a simple differential equation (S) would be valid for /, then also 



pi/' = -E» f {k) 

and taking the square we get 



k=0 



P if ,2 =pi(f-n 



(k) 



J2 Pkf 

k=0 

v mi / 



If we substitute @, and the corresponding representations for (k = 3, . . . , n) which we get by 
differentiating (g), we arrive at an algebraic identity, and the conclusion follows by the transcen- 
dency of /, again. □ 

By similar means we see that the elementary functions cot x, and cscx as well as the corresponding 
hyperbolic functions are not simple. Moreover, the generating functions 

ry. 00 ~k p X/2 °° p h 

of the Bernoulli and Euler numbers Bk, and Ek (k £ No) (see e. g. Q, (23.1)) satisfy the nonlinear 
differential equations 

/'=^((1-^)/-/ 2 ), and f = t-f\ 

respectively, and so by a similar procedure are realized not to be simple. As exactly the simple 
differential equations correspond to simple recurrence equations for the Taylor series coefficients 
(see e. g. @), this implies the following 

Corollary 1 Neither the Bernoulli nor the Euler numbers satisfy a finite homogeneous linear 
recurrence equation with polynomial coefficients. 
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2 Construction of simple functions 

If simple functions are given, then we may construct further examples of simple functions by the 
following procedures: integration, differentiation, addition, multiplication, and the composition 
with certain functions, namely rational functions and rational powers. Before we prove this main 
result, we consider an important particular case. We will prove that if / is simple, and q € N, then 
/ o x 1 ^ is simple, too. To get this result, we need the following 

Lemma 1 Let / : D — > R be infinitely often differentiable, let be D C R + 7^ 0, and let q G 1ST 
with q > 2 be given. Consider the function h : D — > H with h(x) := f(x l l q ). 

(a) Then we have for r G IN with 1 < r < q — 1 

x r/lf(r) {x l/ q) = j2C° qrk X k h^(x) (4) 
k=l 

with constants C® rk , in particular 

Cjn=9, C° qrl =q(q-l)---(q-r + l), and C° rr = <f , 

and 

f^(x^) = j2C 1 g0k x k - 1 h^(x) (5) 
fc=i 

with constants C q0k , in particular 

Cgoi = g! , and = . 



(b) Further for s G IN we have 



/N)( x V 9) = ^q ofe ^-^) (x) 



with constants C q0k , in particular 

C s g ,o, sq = Q sg , (6) 

and for r G IN with 1 < r < g — 1 

sg+r 

x r/qf(.sq+r)( x l/q) = ^ C* rjfc ^ (x) 

fc=s+l 

with constants C| rfc , in particular 

r? s — ,7 s< ?+ r f7i 



4 



Proof: Differentiating the identity /(x 1 ^) = h(x), we get 

x Vif'{xVi) = qx h'{x) , 

a further differentiation yields 

x 2/ q f^h) + = ? 2 x ^ (x) + q 2 x 2 h „ {x) ^ 

and therefore 

x 2/q J"( X V?) = g (g _ 1) a; + g 2 X 2 . 

Differentiating once more and rearranging similarly, gives 

x 3/ q /'"( X V?) = q (g _ i) ( ? _ 2) x + 3 g 2 (g - 1) x 2 /i"(x) + g 3 x 3 ti"{x) . 
Assume now, for r € IN with 1 < r < q — 1 the relation 

k=i 

with 

C° rl = ? (g-l)---( 9 -r + l), and Cj,. = q r 
is valid, then by another differentiation we get 

x (r+l)/q f (r+l)( x l/qj = ( q - r ) C° rl X ti (x) 

= +j2((kq-r) C° qrk + q C^) x k (x) + q C° qrr x r+1 (x) . 

k=2 

For r < q — 1, we set 

Cjr+1,1 == (? - r) C° qrl =q(q-l)...(q- r ), 
Cg,r+l,r+l := 9 Cgrr = Q<f = </* +1 > 

and 

Cj, r+ i >fc := (kq - r) C° qrk + g Cj^^ (2 < A; < r) . 
For r = q the same method applies, if we set 

and 

Cj, fc := ((* - 1) g + 1) Cj,_ 1>fc + qC° m _ hk _ 1 (2<k<q-l), 

which proves (a). 

To prove (b) we assume that for some s G IN 

sq 

f(sq) {x l/q ) = J^ C S okX k-s h (k) {x) 
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with 

C s - n sq 

If we define the functions f s and h s by 

sq 

f s (x) := f&\x) , and h s (x) := ]T C s q0k x k ' s h< k \a 

k=s 

then we may apply (a) , and get for 1 < r < q — 1 that 

x r/lf(r) {x yi ) = j2C° qrk X k hf\x), 
k=l 



and further 



hence 



fM(x 1 /«) = Y,C 1 q0k x k - 1 hW(x), 
k=i 

r 

x r/q f (sq+r) = £ ^ ^fc fe (fc) ( x) > (g) 

fc=l 

/((^)( x V ?) = ^ C i ofe ^-i /l W (x) . (9) 
fc=i 

To finish our proof, we will calculate the derivatives (k = 1, . . . , q) algorithmically. We have 

sq 



and 



h s (x) = J2C s q0j xi- s hd\x) 



J=S 

and therefore 

sq sq 

K(x) = Y, CS q oj^ ShU+1 \ x )+ E (j-s)C s q0j x^ s - 1 h^(x) 

3=8 j=s+l 
sq+1 sq 

= E C s qAj _ 1 xi- s - 1 h«\x)+ J2 (j-s)C s q0j x^ s - 1 h^(x) 

j=s+l j=s+l 
sq 

= E ((i" s ) % + C?,cw-i) 1 ^ (*) + C«,o, S9 (*) • 

j=s+l 

We may write this as 

sq+l 

h' s (x)= D]x j - s - 1 h^\x) 

j=8+l 

with 

,i ._ / U - s ) °qOj + C ? S ,oj-i if s + 1 < 3 < ^ 



D 



Continuing iteratively in the same fashion, we get representations 

sq+k 
j=s+k 



with 



D k l\ + (j - s - k + 1) Z^ -1 ifs + fc<j<sg + fc 



J 1 g S9 if j = s </ + k 

We substitute these equations in (||) to get 

r sq+k 

x r/q f(s q +r)( x l/qj = ^ g Cj* £)* j^"' fcC?) ( x ) . 
fc=l j'=s+fe 

For j with s + 1 < j < s? + r we use the notation 

Mj := {k e INT | 1 < A; < r and s + fc<j<sg + £;} 
Then M S(J+r = {r}, and 

x r/q f (sq+r)( x l/q^ = g CJ^ a^"' fcC?) ( X ) 

with 

qrj ' / y ^qrk j 
keMj 

for s + 1 < j < s q + r, in particular 

g,r,sg+r ^ qrr sq+r H H H 

proving ([?]). 

To prove (||), we substitute Equations (|io|) in @ to get 



/((»+i)9)(a;i/«) = C s +l x j - s ~ l h®(x) 

j=8+l 

with 

/^is+i \ ^ f~il r>k 

U qOj ■- 1^ ^qOk^j 
keMj 

for s + 1 < j < (s + 1) q, in particular 

° ? ,0,(s+l) 9 - U gOg ^sg+g — 9 9 ~ Q 

finishing our proof. 
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Remark 1 We remark that one easily can establish that all constants C qrk are integers. Further, 
we note that using the differential operator 9 q := qx -5- (compare with the proof of ||, Theorem 
8.1) Equation (Q) gets the simple form 

x r/ q f (r) {x i/ q) = 9q y q _ 1){ Q q _ 2 )... (p q _ ( r _ 1)) h ( x ) (1 < r < q ) . 

Next we give some applications of the lemma. First assume that /' = /, a solution of which 
is given by the exponential function f(x) = e x . Then jfw = /, and so f^ q \x l / q ) = /(x 1 ^), and 
therefore we get for h(x) := e xl/q the differential equation 



fc=i 



a/3 



e. g. for g = 3we have the differential equation for e x 

27 x 2 ti" + 54xh" + 6ti = h . 
Further we consider the homogeneous Euler differential equation 

n 
k=0 

with constants a k (k = 0, . . . , n), a n 7^ 0. As by the Lemma we have 

£cj*s fe fc< fc >(x) (l<r<?-l) 



(11) 



t r f (r) (t) 
t q f iq) (t) 
t sq f( sq \t) 



t=x 1 /i 



t=x 1 /i 



t=x 1 /i 



k=l 



k=l 
sq 

J^C s q0k x k h^(x) (s G IN) 



and 



k=s 



sq+r 



t=x 1 /i 

we immediately conclude the following 



C s qrk x k h^{x) (l<r<«f-l) 

k=s+l 



Corollary 2 Assume, the function / satisfies a homogeneous Euler differential equation (O) of 
order n. Then the function / o x l / q satisfies a homogeneous Euler differential equation of the same 
order n. □ 

As an application of Corollary |2| we get 

Corollary 3 Assume P is a polynomial of degree n, then P o satisfies a homogeneous Euler 
differential equation of order < n + 1 . 
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Proof: P satisfies the differential equation x n+l /( ri+1 ) =0. □ 

Our next step is to prove another lemma which we will need to prove that if / is simple, and g£N, 
then so is / o x 1//<? . 

Lemma 2 Let / : D — > H be an infinitely often differentiable function that is simple of order 
n, let be D C H + 7^ 0, and let q £ IN with q > 2 be given. Consider the function h : D — > H 
with h(x) := f(x 1 / q ). Then the linear space Lj over K[x] (i.e. with rational function coefficients) 
generated by the functions x r / q f ( ~ k \x 1 l q ) (r,k <G N ) has dimension <ng. 

Proof: Given simple / of order n, the linear space L/ is generated by the functions of the No x No 
matrix 

f(y), yf(y), y 2 f(y), y r f(y), ■■■ 

f'(y), yf'(y), y 2 f'(y), y r f'(y), ••• 

f"(y), yf"(y), y 2 f"(v), y r /"(y), ••• 

/ (fe) (y), yf {k) (y), y 2 f (k) (y), y r f (k) (y), ••• 

using the abbreviation y := x x l q . 

At most the first q expressions /{x 1 ^), x l l q f{x l / q ), xti-V/i f{x l / q ) of the first row are 
linearly independent (over as the expressions for r = sq (s £ N) depend linearly on the first 

expression f{y) {y sq = x s is a polynomial), further the expressions for r = sq + 1 (s £ N) depend 
linearly on the second expression y /(y), and so on. The same conclusion follows for each other row 
so that we obtain that Lf is generated by the q x No matrix M 

f(y), yf(y), y 2 f(y), y q ~ x f(v), 

f'(y), yf'(y), y 2 f'(y), y^f'iy), 
f"(y), yf"(y), y 2 f"(y), y q ~ l f"(y), 

f (k \y), yf (k) (y), y 2 f {k \y), y"' 1 f (k \y), 



Now, we consider the first column. We will show that at most the first n expressions of this column 
are linearly independent. By the given simple differential equation for / there are polynomials p, 
and Pk (k = 0, . . . , n — 1) such that 



n-1 

p(y)f {n) (y) = Y,PK(y)f (k) (y)> 

k=0 



or 

n-1 

p(x 1 / q )f( n \x 1 / q ) = J2 
k=0 
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turns out to be a polynomial in x, then this equation tells that /^(x 1 ^) depends linearly 
(over -f^[x]) on the elements x T / q f^(x l ' q ) (r £ No, < k < n — 1) of the rows above in our scheme. 
As this, however, is not necessarily the case, we construct a polynomial r such that p(x 1 ^ q ) r(x 1 ^ q ) 
is a polynomial in x. 

Therefore we use the complex factorization p(y) = c(y — y\) (y — y?) - • • (y — y m ) of p, and by 
the identity 

y q -yf = (y-yi)Yl y 3 yf 13 

j=0 

we see that multiplication of p by the polynomial 

m I q— 1 \ 
1=1 \j=0 J 

yields a polynomial s in the variable x q . Even though, for technical reasons, we use a complex 
factorization, it turns out that r £ R[x] (Q[x]) if p € TR[x] (Q[x]). Hence P(x) := s(x 1 ^ q ) = 
p(x l l q )r(x l / q ) is a polynomial in x. Thus we see that multiplying ( |l2l) by r(y) generates the 
representation 

n-l 

P[x) f^ix 1 '*) = p k{* llq ) 

fc=0 

with polynomials P, and Pf- (k = 0, . . . ,rt — 1). 

This equation tells that /( n )(x 1//|J ) depends linearly on the functions of the rows above. Similarly, 
by an induction argument, f^'(x^' g ) (j > n) depend linearly on those functions x r ^ q f( k \x l l q ) 
(0 < r < q — 1,0 < k < n — l)of the first n rows of M. 

Considering the other columns, the same argument can be applied to show that y r /w) (y) depend 
linearly on the functions of the first n rows of M . This shows finally that Lf is generated by the 
nq functions x r l q f <yk \x l l q ) (0 < r < q — 1, < k < n — 1), and therefore has dimension <nq. □ 

Lemma |l] shows immediately that for h := f o x 1 ^ the linear space L/j over K[x] generated by 
h, h! , /i", ... is a subset of Lf (declared in Lemma |2|), and so by Lemma ||| is of dimension < nq. 
Thus we have proved the following Theorem (compare []T5|| , Theorem 2.7, Maple function 
algebraicsubs). 

Theorem 2 Let / be simple of order n, and let Then / o x l / q is simple of order < nq. □ 

Whereas Theorem || gives a complete answer with regard to the existence of a simple differential 
equation, in |jl5|| , Theorem 2.7, it is shown that this is true for the composition with each algebraic 
function (p(x) with (fi(0) =)), it does not tell anything how such a differential equation may be 
obtained. Therefore we continue to describe an algorithm that generates the differential equation 
of lowest order for h, which gives a second proof for Theorem ^. 

Algorithm to Theorem |^: Given the differential equation 

n 

^2pk(x) f {k) {x) = (p k polynomials (k = 0, ... ,n), p n 0) (13) 

k=0 
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for / we have to construct a simple differential equation for h(x) := /(a; 1 ' 9 ) (5 £ K,? > 1). In 
this construction Lemma [l] will play a key role as it did in the proof of Corollary 0. Therefore we 

N 

prepare the given polynomials pk in the following way: If P(x) = a,j x 3 is any polynomial with 

r 

(Z7v 7^ 0, then we decompose P by the q polynomials 

N 

P m (x) := Y a 3 xj (m = 0, 1,... ,g- 1) 

j=0 
j=m (mod q) 

in the form 

q-l 

P(X) = J2 P m{x) , 
m=0 

and P is the zero polynomial if and only if P m is the zero polynomial for all m = 0, 1, . . . , q — 1. If 
m G {0, 1, . . . , q — 1} then j = m (mod q) if and only if j = Sj q + m with Sj € iVo. Therefore we 
have 

N 

Pm(x)=X m Y, ">■'""■ 

j=0 
j=m (mod q) 

If we set 

N 

P* m {y):= Y "if* 

j=0 
j=m (mod q) 

q-l 

we have P m (x) = x m P^ l (x q ), and therefore we arrive at the decomposition P(x) = x m Pm{x q ) 

m=0 

with the polynomials P^. For the degrees of the polynomials Pm(y) we obtain the relations 

!< sn (0 < m < r^) 

= s N ( m = r N ) , (14) 

< sat — 1 (r^v < m < q — 1) 

and P* N (y) ^ 0, where we use the representation N = sjy q+fN (stv £ INoi r 7V £ ]No> < < <Z— 1). 

Now we continue with (|l^). We differentiate (|l~3| ) iteratively Z* := n(q — 1) times to get the 
equations 

YPki(x)f {k) (x) = (15) 

fc=0 

with polynomials (Z = 0, 1, . . . , I*), and we note that p n +i,i = Pn ^ 0, and that p^o = Pk- These 
form a set of I* + 1 = n(</ — 1) + 1 equations. 

We decompose each polynomial pki (Z = 0, 1, . . . , Z*) in the way introduced above to get 

q-l 

Pkl(x) = Y X m Qklm(x q ) 
m=0 
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with polynomials Qklm- This brings our I* + 1 equations (15) into the form 
n+l q— 1 



£ £ Qklm(x q ) (x m f^(x)) =0 (I = 0, 1, . . . , f) . 



(16) 



k=0 m=0 



If we write k = s^q + (sk € XNoj r fc = 0, 1, . . . , q — 1), a brief look in Lemma [l] shows that for 
the purpose to arrive at a differential equation for h the "variables" to be eliminated in the above 
equations are the terms (x m f {k \ x) ) where m G {0, 1, . . . , q — 1} \ {r k }- Every derivative so 
generates q— 1 unknowns, and as we have nq+1 derivatives, these are (nq+l)(q— 1) = q(l* + 1) — 1 
unknowns. Therefore we need to generate equations out of each of our I* + 1 equations ( |l6| ) 
to arrive at a set of g(/* + 1) equations, i. e. enough equations to yield the differential equation 
searched for. 

This is done by multiplying each of the equations (jl^) by the factors x p (p = 0, . . . , q — 1) 
generating the q(l* + 1) equations 



n+l q+p—1 n+l q—1 

E E Qk,i,m- p(* q ) (x m / (fc) (x)) = E£W* 5 ) (W (fe) (z: 

fc=0 m=p fc=0jr'=0 
(I = 0, 1, . . . , T, k = s k q + r k , s k G Ko, r fc = 0, 1, . . . , g - 1, p = 0, 1, 
of the g(Z* + 1) — 1 unknowns f^ k \x)^j , where we use the abbreviations 



0, 



(17) 



Qkipjiy) : ~- 



yQk,i, g +j-p(y) if i = o, ... ,p- i 

Qk,i,j~ P {y) i£j=p,...,q-l 



and 



Qkiojiv) ■= Qkijiy) 



(j = 0,l,...,g-l) . 



(18) 



(19) 



We will show that these unknowns always can be eliminated, therefore arriving at a differential 
equation of order nq for h. 

We note, that if we differentiate (|l3|) iteratively, generating q new equations at each step by 
multiplication with x p (p = 0, . . . ,q — 1), and checking at each step if the unknowns can be 
eliminated, this algorithm obviously results in the differential equation of lowest order valid for h. 

Now we show that in the last step the algorithm always succeeds. Therefore we rewrite the 
equations system (|T7|) using the q x q matrices 



Au{y) ■= (Qkipj 



With the aid of these matrices each block of equations with fixed / can be written as the matrix 
equation 

/ f {k) {x) \ 
xf( k \x) 
(x) 



n+l 
k=0 



x 2 /W(a 



. 



y x q-lf(k)^J 
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Because of ( JTq) and (|T^) the matrices Aki(y) are of the following special form 



A k i(y) 



( Qkio(y) 
yQk,i, q -i(y) 
yQk,i, q -z(y) 

yQm(y) 
yQkn{y) 

n+l 



Qkiiiy) 
Qkio(y) 
yQk,i, q -i(y) 

yQm(y) 
yQkniy) 



Qm(y) 
Qkiiiy) 
Qkio(y) 

y Qku{y) 
yQm{y) 



Qk,l,q-2{y) 

Qk,i, q -z{y) 

Qk,l,q-4(y) 

Qkio(y) 
yQk,i, g -i(y) 



Qk,i, q -i(y) \ 

Qk,l,q-2{y) 
Qk,l,q-3(y) 

Qkiiiy) 
Qkio(y) ) 
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Prom the relations J2 Pm(x) f( k \x) = with p ki (x) = x m J2 Qkim(x q ), p n +i,i(x) = p n (x) 0, 

k=0 m=0 
and Pko(x) = Pk(x) we see that the matrix Aki(y) is completely determined by the polynomial 

Pki(x), and we have 

A n +i,i{y) = A n0 (y) =: A n (y) . 

In particular, the matrix A n (y) is completely determined by the nonvanishing coefficient polynomial 
Pn(x) of f( n ^ in the original equation (|l3|). We will show now that the matrix A n (y) is invertible. 

N 

Let p n {x) = CLjX J with ajy / 0, N = sjy q + (sjv,rjv € No, < < q — 1). We have 

3=0 

q-1 

p n (x) = Pno(x) = x m Qn0m(x q ). If we use the notation Q m (y) := Q n 0m 
degree relations (113) read 



(y), then for Q m the 



deg Q m {y) 



< 



< 



sn (0 < m < rjy) 

sn (m = r N ) 

sn — 1 (rN<m<q — l) 



and we have Q rN ^ 0. This information leads to the following observations about the degrees of 
the polynomial entries of the matrix A n : The entries of A n for which the difference of colomn index 
and row index A is > (r^r + 1) (region I) have degrees < sjv — 1, the entries for which < A < rjy 
(region II) have degrees < sjv, further if — + 1 < A < — 1 (region III), then the degrees again 
are < sn, and finally if A < — (region IV), then the degrees are < sjv + 1. 

Assume next that p n (0) 7^ 0, then we get |^4 n (0)| = Qq(0) = p^(0) 7^ 0, and therefore A n (y) 
is invertible. To obtain the same result if p n (0) = 0, we use our degree observations. By the 
Weierstrafi representation of the determinant D = \djk\ of a q x q matrix (djk) 

D = ^2 si g n ^^1,77(1) ^2,77(2) • • • d qj7r ( q ) 

where the sum is to be taken over all permutations ir € S q , we observe that in the Weierstrafi 
representation of the determinant of our matrix A n (y) the summands 

Q q (y), ±yQ q 1 (y), ±y rN Q q N (y), y^ 1 Q q q _M 

occur. Now it follows that 

deg (y rN Q q rN {y)) =qs N + r N = N = degp n (x) > , 
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and for all other summands P(y) using our degree observations we get the relation 

deg P(y) < qs N + r N = N . 

Hence the degree of |j4 n (y)| equals N, thus |A n (y)| cannot be the zero polynomial, and A n {y) 
therefore is invertible. 

Using the inverse A^ 1 , we can now write 



x /(rH-Q( x ) 



!y(n+Z)( a 



n+l-1 

E A-\x«)-A kl (x«) 

k=0 



( f ik) W \ 
xf( k \x) 

X 



x 2 / (fe) (i 



(z = o,...,r) . 



V x«- 1 /( fe )(»)/ 



Using the fact that the inhomogeneous parts of these equations do not vanish, it is now easily 
established that the unknowns can be eliminated, which finishes the algorithm. □ 



After these preparations we are ready to state our main theorem (compare |L5|], [|14[). 

Theorem 3 Let the two functions / and g be simple of order n and m, respectively, and let r be 
rational. Then 

(a) jf{x)dx is simple of order <n+l, 

(b) /' is simple of order < n , 

( c ) / + 9 is simple of order < n + m , 

(d) f ■ g is simple of order <nm, 

(e) for is simple of order < n , 

(f) / o x p l q (p, q € 2) is simple of order < nq . 

Proof: (a): Let h = J f(x) dx, and let / satisfy the simple differential equation @, then obviously 

n 
k=0 

and so h is simple of order < n + 1. 

(b): Let h = /', and (^) be valid. If po = then (||) is a simple differential equation valid for h, 
and so h is simple of order n. If po ^ 0, however, then, dividing by po, we write (0) as 



fc=i 



(20) 



with rational functions r&. If we differentiate the original differential equation, we get a further 
differential equation 

n+l 

E^-/ ( 

fc=0 



p(fc) 
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with polynomials Pk (k = 0, . . . ,n + 1). We substitute / by pQ), and multiply by the common 
denominator to get a simple differential equation of order < n for h. 

(c): First we give an algebraic argument for the existence statement. Let / and g satisfy simple 
differential equations of order n and m, respectively. We consider the linear space Lf of functions 
with rational coefficients generated by /, /', /", . . . , f^ k \ .... As /,/',... , are linearly depen- 
dent by (|2|), and as the same conclusion follows for /', /", . . . , /( ra+1 ) by differentiation, and so on 
inductively, the dimension of Lf is < n. Similarly L g has dimension < m. We now build the sum 
Lf + L g which is of dimension < n + m. As f+g, (/+ g)', . ■ ■ , (f-\-gp k \ ■ ■ ■ are elements of Lf + L g , 
arbitrary n + m + 1 many of them are linearly dependent. In particular, / + g is simple of order 
<n + m. 

Now we like to present an algorithm which generates the differential equation for the sum: We 
may bring the given differential equations for / and g into the form 



n-l 



/ w = E^/ a) > ( 21 ) 

j=0 



and 

g 



m—l 
(m) _ ST n, 



J2lk9 (k) , (22) 

k=0 

with rational functions pj , and respectively. By iterative differentiation and recursive substitution 
of (pH), and (22), we generate sets of rules 



n-l 

f (l) = J2P l jf (j) (l = n,...,n + m) , (23) 

j'=o 

and 

m—l 

5 (/) =E^5 {fc) (/ = m,...,n + m), (24) 

fc=0 

with rational functions p'-, and q l k , respectively. Now, by iterative differentiation of the defining 
equation h := / + g, we generate the set of equations 

h = f + g 
ti = f' + g' 

h" = f' + g" . (25) 

fo(n+m) _ j(n+m) _|_ g(n+m) 

Next, we take the first max{n, m) of these equations, and use the rules ( |23"|) and ( |24] ) to eliminate 
all occurrences of /W (Z > n), and (Z > m). On the right hand sides remain the n + m variables 

(Z = 0, . . . , n — 1) and (Z = 0, . . . , m — 1). We solve the remaining linear system of equations 
for the variables (1 = 0,..., max{n, m}) trying to eliminate the variables f( l > (I = 0, . . . , n — 1) 
and (Z = 0, . . . , m — 1). If this procedure is successful, it generates the simple differential 
equation for h, searched for. 

If the procedure fails, however, we go on taking the first max{n,m} + 1 of equations (p5|), 



doing the same manipulations, and so on, until we take the whole set of equations (25), where the 
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procedure must stop. 

(d) : This situation is handled in exactly the same way as case (c). The only difference is that to 
differentiate the product h := fg we use the Leibniz rule, and the variables to be eliminated are 
the n m products /w) g( k ) (J = 0, . . . , n — 1, k = 0, . . . , m — 1). It is easily seen that this procedure 
generates a differential equation for h of order <nm. 

That the algorithm stops is seen by the fact that finally we arrive at a set of n m + 1 equations 
with nm variables to be eliminated, where all equations actually possess a nonvanishing inhomo- 
geneous part 

(e) : Let h := f o r for some rational r, and let / satisfy dH). We compose (g) with r, and get 

n n 

= ^2{pk°r) • (/^or) = X! * ' (f^ or ) ilk rational functions) . (26) 

k=0 k=0 

Differentiating the identity h = f or leads to 

h = f or • r or f or 
h" = for ■ r" + for ■ r' 2 or for 

and inductively it follows 

k 

^(k) _ . |j(i)orj (Vj rational functions) . 

i=i 

Solving for f^or and substituting the previous results we are lead to a representation 

fWor = y £R J -h^ 
3=1 

with rational functions Rj. 

Substituting these results in ( f26| ) and multiplying by the common denominator yields a simple 
differential equation for h, and we see that the order of h is < n. 

(f) : The statement follows immediately by an application of Theorem || and (e). □ 

Remark 2 We remark that again our proofs provide algorithms to generate the differential equa- 
tions searched for. 

Further we note that by Theorem [l] the functions 1/g and f/g in general are not simple, if / 
and g are, as secx = 1/ cosx and tanx = sinx/ cosx show. 
Note that Theorem || (f) generalizes Theorem 8.1 in ||. 

The following consideration strengthens Theorem ||| (c) , and brings it in connection with the fun- 
damental systems of simple differential equations. 

Corollary 4 Assume, / satisfies a simple differential equation F of order n, g satisfies a simple 
differential equation G of order m, and F and G do not have any common nontrivial solution. 
Then the order of any simple differential equation satisfied by the sum h := / + g is > n + m. In 
particular, h satisfies no simple differential equation of order < n + m, and exactly one differential 
equation H of order n + m which is generated by the algorithm given in Theorem ^ (c) . 
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Proof: Let fx , fa, . . . , f n be a fundamental system of the differential equation F, and gi , , . . . , g m 
be a fundamental system of the differential equation G. Then fx, f2, ■ ■ ■ , fn are linearly independent 
(in the usual sense), and gi, 52, • • • , 9m are linearly independent, as well. Assume now, the n + m 
functions fx, /2, ■ ■ ■ , f n , 91,92, • • • ,9m are linearly dependent. Then there are constants dj (k = 
j, . . . ,n) and bk [k = 1, . . . , m), not all equal to zero, such that 

n m 

E«i/i + E 6 ^fc = 0. (27) 
3=1 k=i 

From the linear independence of the subsets fx, /2, . . . , f n , and g±, g2, ■ ■ ■ , g m it follows that at least 
one of the numbers 

a, (j = 1, . . . , n), and one of the numbers bk (k = 1, . . . , m) does not vanish. Therefore by the 

n m 

linearity of F and G the functions f '■= a j fj and g '■= — bj~ gk are solutions of F, and G, 

3=1 ' fc=i 

respectively, which by (|27| ) agree, and we have a contradiction. 

Thus we have proved that /1, /2, • • • , f n , 9\,92, ■ ■ ■ ,9m are linearly independent. As the zero 

function is a solution of both F, and G, it turns out that fi, f2, ■ ■ ■ , fn, 9i, 92, ■ ■ ■ , 9m must satisfy 

any linear differential equation H for the sum, and so form a fundamental system for H. Therefore 

the order of H is > n + m. □ 

We mention that, on the other hand, if the algorithm presented in Theorem [|| (c) generates a 
differential equation of order < m + n for h := / + g, then the simple differential equations F for 
/ and G for g must have a common nonzero solution, i. e. the fundamental systems F and G are 
linearly dependent. We gain this insight without solving any of the differential equations! 

From Theorem |3| (d) a similar statement can be obtained for the product h := fg: If our 
algorithm generates a differential equation of order < nm for h, then the product of the fundamental 
systems Sf and Sg of F and G, i. e. the set 

{f 9 1/ £ Sf, g £ <Sg}, is linearly dependent. Again, we gain this insight without solving any 
of the differential equations. 

In connection with Theorem ||| (d) for the product h := fg we note that for the special case of 
the square h := f 2 = f ■ f the bound for the order of h can be considerably strengthened. 

Corollary 5 Let / be simple of order n, then h := f 2 is simple of order < 5^±il, 

Proof: A careful study of the algorithm given in the proof of Theorem ^ (d) shows that the order 
cannot exceed niZLtil. This depends on the fact that the variables to be eliminated are the n ^" 2 +1 ^ 
different products ft) (J = 0, . . . , n - 1, k = 0, . . . , n - 1). □ 



3 Algorithmic calculation of simple differential equations 

Theorem |3| enables us to define the linear space F of functions generated by the algebraic functions, 
exp x, In x, sin x, cos x, arcsin x, arctan x, (and, if we wish, further special functions satisfying simple 
differential equations,) and the functions that are constructed by an application of finitely many 
of the procedures (a)— (f) of Theorem ||. The theorem then states in particular that F forms a 
differential ring. 

The algorithm which generates the simple differential equation of lowest order for / given in 
H U is strengthened by Theorem S 
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Algorithm 1 (Find a simple differential equation) Let f £ F. Then the following procedure 
generates the simple differential equation of lowest order valid for /: 

(a) Find out whether there exists a simple differential equation for / of order N := 1. Therefore 
differentiate /, and solve the linear equation 

f'(x) + A f(x) = 

for Aq; i. e. set Aq := — y^y- Is Aq rational in x, then you are done after multiplication with 
its denominator. 

(b) Increase the order N of the differential equation searched for by one. Expand the expression 

/W(s) + A^f (x) + ■■■ + A f(x) , 

and check, if the summands contain exactly N rationally independent expressions (i. e.: 
linearly independent over iif[z]) considering the numbers Aq, Ai,... , Ajy-i as constants. Just 
in that case there exists a solution as follows: Sort with respect to the rationally independent 
terms and create a system of linear equations by setting their coefficients to zero. Solve this 
system for the numbers Aq, Ai, . . . ,An—i- Those are rational functions in x, and there exists 
a unique solution. After multiplication by the common denominator of Aq, Ai, . . . , An-i you 
get the differential equation searched for. Finally cancel common factors of the polynomial 
coefficients. 

(c) If part (b) was not successful, repeat step (b). 

Proof: The proof given in || shows that the given algorithm indeed generates the simple differ- 
ential equation of lowest order valid for / whenever such a differential equation exists. Theorem || 
now guarantees the existence of such a differential equation, and therefore the algorithm does not 
end in an infinite loop. □ 



Remark 3 Obviously / € F can be checked by a pattern matching mechanism applied to the 
given expression /, and Theorem || then moreover gives a priori estimates for the possible order of 
the solution differential equation. 

We further point out that an actual implementation of the algorithm should be able to decide 
the rational dependency of a set of functions. Otherwise, it may fail to find the simple differential 
equation of lowest order. An example of this type is given by f(x) := sin (2x) — 2 sin x cos x, for 
which the algorithm yields the differential equation f" + / = rather than /' = if the rational 
dependency of sin (2a;), and sin x cos x is not discovered. 

In a forthcoming paper [jll]] we discuss the more general situation of functions / that depend on 
other special functions not contained in F, like Airy functions, Bessel functions and orthogonal 
polynomials. 

We mention that the algorithms presented in Theorem ||| can be combined to get the following 
different 

Algorithm 2 (Find a simple differential equation) Let be / £ F. Then recursive application 
of the algorithms presented in Theorem ^ generates a simple differential equation for /. 
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Remark 4 We note that, in general, we cannot control the order of the resulting differential 
equation as the algorithm uses recursive descent through the expression tree of the given /. Given 
a sum of two functions or order n and m, however, it is easily seen that the order of the resulting 
differential equation is the lowest possible order which is > max{n, m}. So with the given algorithm, 
it is principally impossible to derive any differential equation of order lower than two for the example 
function f(x) := sin (2x) — 2 sin x cos x as the lowest order differential equations for the summand 
functions sin (2x), and sin x cos x are of order two. 

This can be interpreted as follows: By construction, for a sum / := /1 + /2 a differential equation 
is obtained that is valid not only for f\ + f 2 , but for any linear combination Ai f\ + A2 ji or j\ 
and ji. In other words, this subalgorithm generates a differential equation for the linear hull of f\ 
and f 2 . 



through the quotient 



4 Sharpness of the orders 

In this section we give examples that show that the bounds for the least orders that are given in 
Theorems [2]-|3], and in Corollary [| in fact may be assumed. 
First we consider the statement of Corollary |5[ We get 



Theorem 4 For each n £ N there is a function f n that is simple of order n, such that h := is 

n+ 
2 



simple of order MlLtl) Indeed, an example function of that type is 



" h 



/„(*):= 5> (28) 



k=l 



Proof: We show that f n , given by (|28|), is simple of order n, and that f% is simple of order "( n+1 ) 



2 

k 

As each function e x [k = 1, . . . , n) is simple of order 1 (satisfying the differential equation 
/' — kx k ~ 1 f = 0), by an inductive application of Theorem || (c), we see that f n is simple of order 
< n. To show that f n is simple of order > n, we first show that e x (k = 1, . . . ,n) are rationally 
independent. Suppose, any linear combination 



ne x + r 2 e x +--- + r n e x =0 (29) 



with rational functions (k = 1, ... ,n) representing zero is given. From ( |29|) we get, assuming 



r k e x =- E 



n e 



k 

and dividing by e x , we get 



e 



19 



with rational functions Rj (j = 1, . . . , n). For k = n, we get by taking the limit x — > oo, along the 
positive real axis 



n-l 



1 = - lim y R je x '- X " = , 

a;— >oo ' J 

a contradiction to the assumption r n ^ 0. For k < n, the right hand limit 



lim V , ' ' 



does not exist, whereas the left hand side equals 1, a contradiction to the assumption r& ^ 0. 
Therefore rjt = (k = 1, . . . , n), and e x (k = 1, . . . , n) are rationally independent. 

From this it follows that applying Algorithm [l] to f n , in each step n rationally independent are 
involved (the derivative of e x is of the same type), so that the algorithm can be successful not 
earlier than in the nth step. As this algorithm always finds the simple differential equaton of least 
order, the order of f n is > n. 

Next we consider the function Expanding f% we get the representation 



X" 
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/2 (x) = e 2* + e 2x< + . . . + e 2x" + 2 e x f + . . . + 2 e x"- e 

consisting of n+ ^ ^ ) = w ^" 2 f1 ^ summands. An application of Theorem ||| (c), again, shows that /, 
is simple of order < 2" , as each summand / := e xi e xk satisfies the simple differential equation 

f'-2(jxi + kx k )f = 

of order 1. Similarly as above one realizes that e xj e xk (j,k = 1, ...,ra (j 7^ k)) are rationally 

independent, and therefore by an application of Algorithm [l] we conclude that is simple of order 

> Ei!Lhil ) finishing the proof. □ 

Next we consider the statement of Theorem |2[ We get 

Theorem 5 For each n, q £ N there is a function f nq that is simple of order n, such that h := 
/ o x x l q is simple of order n q. Indeed, an example function of that type is 

n 

f nq (x):=Y,e xkq+1 ■ (30) 

k=l 

Proof: Everything follows as in Theorem 0, if we show that 



/( x ) :=e ^ +1 {k = l,...,n, q € N) 



is simple of order 1, whereas 

h(x) := e xk+1/q (jfe = l,...,n, q€N) 

is simple of order q. The first statement was already proven so that it remains to prove the second 
statement. This, however, turns out to be a consequence, of Lemma |l] (a). We start with one 
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rationally independent function, namely h itself. As h(x) = f{x 1 ^ q ), Lemma [l] (a) shows that in 
the rth differentiation step (r < q) for h the new summand 

x r/qf(r)( x l/q) 

is introduced. By a similar argument as in Theorem |] one sees that for each r < q — 1 the functions 
involved functions 

[h(x), x 1/q f(x 1/q ), . . . , x r l q / (r) (x 1/9 )} 

are rationally independent. 

In the qth step, however, by (|5|), rationally depends on h,h' , . . . , M 9-1 **, so that an arbitrary 
linear combination of h,h' , . . . contains exactly q rationally independent terms, showing that h 
is of order q. □ 

Finally we consider the statements of Theorem ||. We get 

Theorem 6 The statements of Theorem |3| are sharp which is seen by the following example func- 
tions: 

(a) For each n the function f n given by (p§|) is simple of order n, and / f n is simple of order n+1. 

(b) For each n the function f n given by fl2~§D is simple of order n, and f' n is simple of order n, too. 

(c) For each n, m € IN there are functions / and g that are simple of order n, and m, respectively, 
such that h := / + g is simple of order n + m. Indeed, example functions of that type are 

n n+m 

f(x):=J2e xk and g(x) := £ ^ • ( 31 ) 

fc=l fc=n+l 

(d) For each n, m € IN there are functions / and 5 that are simple of order n, and m, respectively, 
such that h := / g is simple of order nm. Indeed, example functions of that type are given 

by (1). 

(e) For each n the function f n given by (p£|) is simple of order n, and for each rational function 
r the function / o r is simple of order n, too. 



(f) For each n the function f nq given by (|3y) is simple of order n, and for each p £ IN the function 
/ o x p/i is simple of order nq. 

Proof: (a) By Theorem |] f n is simple of order n. As the antiderivative of f n is not an elementary 
function, it is rationally independent, and so its order is >n+l. 

(b) As the derivatives of each summand of f n depend rationally on f n , by the linearity the order 
of f' n equals n, too. 

(c) In Theorem |] we already proved that / is simple of order n, and that / + g is simple of order 
n + m. A similar argument shows that g is simple of order m, and we are done. 

(d) If we expand the function / g, we get 

n n+m n m 

f(x)g(x) = Y,e xk £ e * k = STJ2e xk e xJ+n 

k=l k=n+l k=l j=l 

with nm rationally independent expressions. Thus the order of / g is >nm. 

(e) By the argumentation of Theorem one sees that f n o r has order n. 

(f) This follows from (e) and Theorem 0. □ 



21 



References 

[I] Abramowitz, M. and Stegun, I. A.: Handbook of Mathematical Functions. Dover Publ., New 
York, 1964. 

[2] Chudnovsky, D. V. and Chudnovsky, G. V.: On expansion of algebraic functions in power and 
Puiseux series I. Journal of Complexity 2 (1986), 271 - 294. 

[3] Chudnovsky, D. V. and Chudnovsky, G. V.: On expansion of algebraic functions in power and 
Puiseux series II. Journal of Complexity 3 (1987), 1 - 25. and products. 

[4] Gruntz, D. and Koepf, W.: Formal Power Series. Konrad-Zuse-Zentrum Berlin (ZIB), Preprint 
SC 93-31, 1993. and 

[5] Hearn, A.: Reduce User's manual, Version 3.4. The RAND Corp., Santa Monica, CA, 1987. 

[6] Koepf, W.: Power series in Computer Algebra. J. Symb. Comp. 13 (1992), 581-603. 

[7] Koepf, W.: Algorithmic development of power series. In: Artificial intelligence and symbolic 
mathematical computing, ed. von J. Calmet and J. A. Campbell, International Conference 
AISMC-1, Karlsruhe, Germany, August 1992, Proceedings, Lecture Notes in Computer Science 
737, Springer- Verlag, Berlin-Heidelberg, 1993, 195-213. 

[8] Koepf, W.: Examples for the algorithmic calculation of formal Puiseux, Laurent and power 
series. SIGSAM Bulletin 27 (1993), 20-32. 

[9] Koepf, W.: A package on formal power series, Konrad-Zuse-Zentrum Berlin (ZIB), Preprint SC 
93-27, 1993. 

[10] Koepf, W.: A new algorithm for the development of algebraic functions in Puiseux series. 
Preprint B-92-21 Fachbereich Mathematik der Freien Universitat Berlin, 1992. 

[II] Koepf, W.: Algorithmic work with orthogonal polynomials and special functions. In prepara- 
tion. 

[12] Maple: Reference Manual, fifth edition. Watcom publications, Waterloo, 1988. 

[13] Neun, W.: Reduce package FPS. Konrad-Zuse-Zentrum fur Informationstechnik Berlin, 1994, 
in preparation. 

[14] Salvy, B. and Zimmermann, P.: GFUN: A package for the manipulation of generating and 
holonomic functions in one variable. Rapports Techniques 143, INRIA, Rocquencourt, 1992. 

[15] Stanley, R. P.: Differentiably finite power series. Europ. J. Combinatorics 1, 1993, 175-188. 

[16] Szego, G.: Orthogonal Polynomials. Amer. Math. Soc. Coll. Publ. Vol. 23, New York City, 
1939. 

[17] Tricomi, F. G.: Vorlesungen uber Orthogonalreihen. Grundlehren der Mathematischen Wis- 
senschaften 76, Springer- Verlag, Berlin-Gottingen-Heidelberg, 1955. 



22 



[18] Wall, B.: Puiseux expansion. An annoted bibliography. Technical report, RISC-Linz Report 
Series No. 91-46.1. 

[19] Wolfram, St.: Mathematica. A system for doing mathematics by Computer. Addison- Wesley 
Publ. Comp., Redwood City. 



23 



